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结构 的 追踪 模型 


间 纵 向 关系 具有 重要 作用 ， 也 为 因果 关 


多 式 的 模型 ， 如 何 


选择 


适当 的 模型 是 重要 的 议题 。 本 文 对 各 模型 进行 概述 ， 并 从 模型 结构 、 预 设 轨迹 、 时 间 点 要 求 
等 方面 进行 比较 ， 最 后 通过 一 个 实例 说 明 如 何 选择 适当 的 模型 。 结 果 表 明 , 不 同 模型 在 变量 
关系 的 判断 上 可 能 给 出 很 不 同 的 结果 ， 实 际 运用 中 应 当 有 模型 选择 和 模型 比较 的 意识 。 
关键 词 ”纵向 关系 交叉 滞后 追踪 模型 模型 选择 
1 引言 

追踪 设计 的 广泛 使 用 推动 了 心理 学 等 社会 科学 领域 对 变量 之 间 因 果 关 系 的 讨论 , 有 学 者 


认为 基于 追踪 数据 建立 的 时 序 关 系 本 质 上 是 相关 性 的 ， 不 足以 支撑 可 靠 的 因果 推论 (Usami 
etal., 2019)， 因 果 推 论 还 是 需要 有 理论 、 文 献 或 经 验 常 识 的 支持 ( 温 忠 记 ,2017)。 因 此 ， 研 


完 者 们 在 使 用 追踪 数据 探究 变量 


% (reciprocal relation) 而 非 
von Eye, 2020)。 纵 向 相互 关系 是 一 种 双向 的 效应 ， 在 实证 
响 另 一 个 变量 的 后 期 取 值 〈 或 变化 量 


Wiedermann & von Eye, 2020)。 这 种 相互 作用 在 心理 学 研究 中 很 常见 ， 例 如 发 展 心理 入 
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间 关 系 时 更 多 使 用 纵向 关系 (longitudinal relation) 或 相互 关 


因果 关系 (causality) 的 表述 (Usami etal., 2019; Wiedermann & 


完 中 


HAE 


{d 


表现 为 一 个 变量 的 


a), PEZ INP (Curran et al., 


2014; 


学 的 理 


论 通常 认为 个 体 特征 变量 与 环境 变量 之 间 是 相互 影响 的 《Curran etal.,2014)。 探 究 变量 间 的 
纵向 相互 关系 对 于 理解 生命 历程 中 的 个 体 发 展 非常 重要 ， 也 能 为 因果 推论 的 验证 奠定 基础 
”本 研究 得 到 国家 自然 科学 基金 项 目 (31771245) 资 助 。 
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(Mund & Nestler, 2019). 


交叉 滞后 面板 模型 (Cross-Lagged Panel Model, CLPM) 是 公认 的 探究 变量 之 间 动 态 关 系 


的 强 有 力 的 方法 。 与 传统 的 分 析 重 复 测 量 数据 的 方法 相 比 ，CLPM 具有 高 度 的 灵活 性 ， 能 够 


BE 


处 理 时 变 协 变量 、 多 变量 分 析 以 及 多 组 分 析 等 ， 且 可 以 提供 多 种 模型 拟 合 指数 和 修正 指数 


(Curran etal.,2010)。CLPM 通过 交叉 滞后 路 径 反 映 变 量 间 的 历时 性 作用 ， 即 构建 一 个 变量 
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的 先前 水 平 对 该 变量 当前 水 平 的 作用 路 径 〈 称 为 自 回归 效应 ) 以 及 其 对 另 一 变量 当前 水 平 的 


作用 路 径 《〈 称 为 交叉 滞后 效应 )。 交 叉 清 后 结构 的 一 个 突出 优点 是 能 够 控 表 


结果 变量 的 先前 


水 平 ， 即 控制 了 一 种 潜在 的 “第 三 变量 ”(third variable) 的 影响 (Cole & Maxwell, 2003). 


mda Hl] Hc 影响 。 


在 纵向 研究 中 , 结果 变量 的 先前 水 3 


E 


主 同时 和 结果 变量 当前 取 值 以 及 预测 变量 相关 , 因此 


D. 


基于 交叉 滞后 结构 出 现 了 一 系列 追踪 模型 ， 如 随机 截 距 交 叉 清 后 面板 模型 Random 


Intercept Cross-Lagged Panel Model, RI-CLPM)、 结 构 化 残 差 潜在 曲线 模型 (Latent Curve Model 


with Structured Residuals, LCM-SR) 和 潜在 变化 分 数 模型 (Latent Change Score Model, LCS) 


AR. 这些 模型 被 广泛 地 运用 于 变量 间 纵 向 关系 的 探究 。 不 同 模型 各 有 所 长 ， 所 适用 的 数据 结 


构 也 不 同 ， 没 有 任何 一 个 模型 是 万 能 的 (Curran et al., 2014; Mund & Nestler, 2019)。 然 而 ， 


zn 


ir 多 


究 者 并 不 了 解 各 模型 之 闻 的 差异 ， 在 实际 运用 中 也 缺乏 模型 选择 和 模型 比较 的 意识 。 
本 文 对 基于 交叉 灌 后 结构 的 追踪 模型 进行 分 析 和 比较 , 首先 对 上 述 模 型 进行 概述 , UGG 


yx 
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ae 


各 个 模型 之 间 的 差异 和 关联 ， 并 明确 不 同 模 型 所 适用 的 研究 情境 , 最 后 通过 一 个 实例 说 明 如 


何 选择 适当 的 模型 。 


2 模型 概述 
2.1 交叉 滞后 面板 模型 CCLPMD 


Jóreskog (1970) 将 时 间 序 列 分 析 中 的 向 量 自 回 归 模 型 纳入 结构 方程 建 模 框 架 ， 从 而 获 


得 了 CLPM， 其 假设 变量 的 观测 值 


对 于 变量 x 和 y， 被 试 i Gels 


m 


围 


其 中 ，jxe、Jyt 表 示 组 均值 ，xii、 潜 表示 时 间 特 定 的 仿 移 量 。 首 先进 行 中心 化 处 理 ， 提 取出 
时 间 特 定 的 偏 移 量 ， 以 表示 每 个 被 试 在 变量 上 随时 间 的 变动 ， 


绕组 均值 上 下 波动 。 


ND 在 时 间 点 上 (CE1…… TO. 的 观测 值 分 别 为 xjt 和 yt: 
Xit = Hyt + Xi 


yit = Hyt + Yit 


* ”一 
Xit = Xit — Hxt 


Vie = Vit — Hyt 


然后 构建 交叉 滞后 方程 : 
Xit = BaXig-i) T YxtYi(-1) + dxit m 
Vit = BytYiqe-1) + YytXig-i) + dyit 
其 中 ，Bit、Byt 为 自 回归 系数 ，yxe、yyt 为 交叉 滞后 系数 ，dxit、dyit 为 回归 残 差 。 一 般 所 称 
的 交叉 滞后 结构 默认 同时 包含 自 回 归 路 径 和 交叉 滞后 路 径 。 
自 回 归 系 数 和 交叉 滞 后 系数 都 反映 了 被 试 内 效应 〈Curran etal.,2014)， 其 中 ， 自 回归 系 
数 (Bees Bye) 体现 了 变量 的 等 级 稳定 性 Crank-order stability)， 即 被 试 们 在 给 定 变量 上 的 高 
低 顺序 不 随时 间 变 化 的 性 质 ， 交 叉 滞后 系数 Vee Vye) 表示 在 控制 了 上 一 时 刻 的 结果 变量 
bn D 对 自身 的 影响 后 ， 预 测 变量 〈 如 xie_D ) MARAE Cülyno 的 历时 性 影响 ， 
滞后 系数 通常 被 认为 是 推断 变量 间 纵 向 关系 的 关键 所 在 。 
CLPM 假设 不 存在 稳定 的 被 试 间 差 异 , 参数 中 同时 混合 了 被 试 内 效应 和 被 试 间 效 应 ,这 
可 能 导致 错误 的 结论 (Mund & Nestler, 2019). 
2.2 随机 截 距 交叉 滞后 面板 模型 CRI-CLPM) 
Hamaker ÆA (2015) Æ CLPM 中 引入 了 随机 截 距 这 一 潜 变 量 ， 得 到 RI-CLPM。 该 模 


型 假设 变量 的 观测 值 围 绕 自 身 的 特质 性 水 平 上 下 波动 。 


Xit = Hyt + lyi + Xj (4) 
Yit = Hyt + lyi + Vit 
与 公式 COD 相 比 ， 公 式 (4) 增加 了 特质 因子 I、Ly:， 表 示 稳 定 的 特质 性 因素 ， 取 值 随 个 体 
变化 。 特 质 因子 的 所 有 负荷 都 固定 为 1， 也 称 为 随机 截 距 项 (Hamaker etal., 2015 )。 在 分 离 
出 组 均值 和 稳定 特质 因子 之 后 ，RI-CLPM 同样 采用 时 间 特 定 的 偏 移 量 Gens yn? 来 构建 交 
叉 滞 后 模型 ( 同 公 式 3)。 相 比 于 CLPM, RI-CLPM 更 能 准确 地 反映 变量 的 被 试 内 变化 过 程 。 
2.3 结构 化 残 差 潜在 曲线 模型 CLCM-SR) 
2.3.1 潜在 曲线 模型 (LCM) 
结构 化 残 差 潜在 曲线 模型 CLCM-SRO 是 以 潜在 曲线 模型 (Latent Curve Model, LCM) 


为 基础 建构 的 。LCM 也 称 为 潜在 增长 模型 (Latent Growth Model, LGM) 和 潜在 增长 曲线 模 


#4 (Latent Growth Curve Model, LGCM )。 该 模型 由 Meredith 和 Tisak (1990) 提出 ， 其 使 用 
特定 时 间 点 的 观测 变量 来 推断 变量 的 洪 在 增长 趋势 。 


Xit = lyi 十 ArSxi 十 Xi 


(5) 
Vie = lyi + AtSyi + Vit 


Lá = Mor + Gai 
Iyi = Hyi + by (6) 
Sxi = Hys + Çxsi 
Syi = Hys + ysi 
FH, Lg ER x、y 的 初 值 〈 截 距 )，Sw 、Syi 表 示 变 化 率 ，4 表 示 时 间 的 函数 值 〈 例 如 
对 于 线性 增长 模型 ，1 = t 一 1，t > 2); 六 、 志 表示 残 差 ， 同 一 个 变量 的 各 个 时 间 点 的 残 
差 相 互 独立 ， 不 同 变量 的 残 差 在 同一 时 间 点 有 相关 ; Hur、Jx 表 示 所 有 被 试 的 初 值 的 均值 ; 
Hass hys RARAWA RERE eri yio Cesi Gysi EIRE. 
LCM 的 主要 目的 在 于 获得 变量 发 展 轨迹 的 特征 ， 初 值 和 变化 率 可 统称 为 增长 因子 
(growth factors)， 然 而 增长 因子 无 法 反映 两 个 变量 随时 间 协 同 发 展 这 一 被 试 内 过 程 ， 因 此 
这 里 只 考虑 结构 化 残 差 潜在 曲线 模型 (CLCM-SRO. 
2.3.2 结构 化 残 差 潜在 曲线 模型 CLCM-SR) 
Curran S A. (2014). 提出 的 LCM-SR， 是 在 LCM 基础 上 引入 交叉 滞后 结构 ， 使 用 分 离 
出 增长 因子 后 的 残 差 来 建立 交叉 滞后 方程 。LCM-SR 假设 变量 x 和 y 具有 系统 的 随时 间 变 化 
的 趋势 , 并 通过 增长 因子 Gui, Iyi Sy, Sy 反映 。 首 先 将 公式 (5) 中 的 残 差 移 到 方程 左 侧 : 


No 
oot Xit = Xit 一 li = AtSyi 


Ql (7) 
> yit = Yit 一 Lyi — AtSyi 
< 然后 再 构建 交叉 滞后 方程 ， 从 而 实现 残 差 的 结构 化 ， 

c Xit = BaXig-i) T YxtYi(t-1) + dxit 


-一 (8) 
F yit = BytYi-1) 十 YytXi(t-1) + dyit 


在 LCM-SR 中 ， 描 述 变 量 发 展 轨 迹 特征 的 增长 因子 并 不 参与 交叉 滞后 方程 ， 后 者 中 的 
关键 参数 仅仅 是 针对 有 具体 时 间 点 的 ， 这 样 就 清晰 地 分 离 了 被 试 间 差异 和 被 试 内 效应 《Curran 
et al, 2014)。 然 而 ， 有 学 者 认为 LCM-SR 在 构建 纵向 关系 方程 之 前 对 变量 进行 了 去 趋势 处 
FR, 排除 了 增长 因子 的 作用 , 因此 除非 有 足够 理由 认为 变量 的 发 展 轨迹 主要 是 由 时 变 的 遗漏 
变量 导致 的 ， 否 则 使 用 该 模型 可 能 造成 估计 偏差 〈Usami etal., 2019). 

2.4 潜在 变化 分 数 模型 (LCS) 


潜在 变化 分 数 模型 ， 也 称 为 潜在 差异 分 数 模型 (Latent Difference Score Model, LOSM), 


HH McArdle 和 Hamagami (2001) 提出 ，LCS 基于 变化 量 (change) 进行 建 模 ， 其 重点 关注 
一 个 变量 的 先前 水 平 对 于 另 一 个 变量 的 变化 量 的 影响 。 
和 上 述 模型 不 同 ，LCS 考虑 了 测量 误差 ， 
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Xit = futt Exit 


(9) 
Yit = fyit + Eyit 
其 中 , frit» fit 表示 潜在 分 数 ， Exito 2 表示 测量 误差 。 
可 获得 潜在 变化 分 数 ， 
Afxit = fxit — fxi-1) 
(10) 


Afyit = fyit J fyit-1) 
其 中 ，4fit、4fyit 表 示 潜 在 变化 分 数 ， 即 t — 1 时 间 点 和 t 时 间 点 之 间 潜 变量 的 变化 量 。 最 后 
基于 潜在 变化 分 数 构 建交 叉 清 后 方程 ， 


Afrit = Gxi + Bat frice—1y + Yxtfyitt-1) + dxit 


(11) 
Afyit = Gyi + Bytfyice—1) + Yytfxi(t-1) + dyit 


其 中 , Gxi 和 Gyi; 描 述 了 时 间 历 程 中 变量 的 总 体 增长 或 减少 , 称 为 恒定 变化 (constant change), 
而 自 回 归 和 交叉 滞后 部 分 则 为 比例 变化 Cproportional change). LCS 模型 可 以 同时 探究 变量 
长 期 性 的 变化 和 即时 性 的 波动 。 
述 各 个 模型 在 参数 估计 过 程 中 均 默 认 要 求 测量 间隔 等 距 ， 以 保证 有 意义 的 参数 解释 
(Mund & Nestler, 2019). 
3 模型 比较 
基于 交叉 清 后 结构 的 各 个 追踪 模型 结构 相似 且 存 在 一 定 关 联 ( 见 图 1 )。 


(SRA. TERRI TO 


交 

变 式 Y 
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| ^| LCM-SR 到 


(ZEN AGAR. WERT) 


_ 结合 ~ 


Dn mag 


《增长 国 子 》 
图 1 各 模型 间 关 联 


首先 ， 对 CLPM 进行 一 定 的 变换 ， 可 获得 其 他 形式 的 交叉 滞后 模型 。 以 CLPM 作为 基 
准 模型 ， 若 在 方程 中 纳入 取 值 随 个 体 变 化 的 特质 因子 (Ii, 1:2 可 得 到 RI-CLPM: FRH 
LCM 模型 的 方式 分 离 出 增长 因子 (Gui, Ly; Sxi，Syi)， 再 基于 残 差 构建 交叉 滞后 路 径 可 得 到 
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LCM-SR; 若 对 观测 值 不 进行 中 心 化 处 理 而 是 考虑 测量 误差 《如 公式 9)， 然 后 ， 采 用 潜在 分 
数 的 差 值 构建 交叉 消 后 路 径 ， 并 加 入 取 值 随 个 体 变化 的 恒定 变化 因子 (Gx:、Gyi)， 可 得 到 


LCS (Usami et al., 2015)。 因 此 ， 可 以 认为 ，LCM-SR、RI-CLPM fil LCS 都 是 CLPM 在 一 
定 变化 条 件 下 的 特殊 形式 。 

上 述 模型 之 间 存 在 数理 关联 。Usami “A (2019) 归纳 了 一 个 统一 的 结构 框架 ， 该 框架 
由 三 种 通用 方程 构成 ， 分 别 是 测量 方程 、 分 解 方程 和 动态 方程 。 假 设 自 回归 系数 和 交叉 滞后 
系数 均 不 随时 间 变 化 ， 则 三 种 方程 分 别 为 : 


Xit = frit + Exit? Vie = fyit + Eyit (12) 

fiit = Mar + lxi + (t = DS + filer fyie = Wye thi + (= DSi + fe (13) 

fis = Gri + Befa Yahaan + Unie fre = Gr + By piety Fes + ye (04) 

需要 注意 ,为 了 能 够 全 面 地 涵盖 各 个 模型 ， 在 测量 模型 中 考虑 了 测量 误差 ， 引 入 了 潜在 

分 数 (far)， 因 此 ， 为 了 在 随后 的 分 解 方程 中 保持 左右 两 侧 符号 的 一 致 性 ， 采 用 谨 .、 所 .来 
表示 时 间 特 定 的 偏 移 量 或 残 差 ， 对 应 于 1) 那样 的 方程 中 的 xx Yio 

通过 在 三 种 通用 方程 中 删除 和 保留 特定 元 素 ,可 以 获得 不 同 的 交叉 滞后 模型 。 例 如 , 在 

量 方程 中 删除 测量 误差 (eu ), 在 分 解 方程 中 删除 截 距 和 变化 率 因子 ([ru + (t 一 DS. 

同时 在 动态 方程 中 仅 保留 交叉 滞后 结构 与 残 差 CB Kieay  Yefgiac + deel» TRE 


P 


Jl 


CLPM; 基于 CLPM， 在 分 解 方程 中 保留 截 距 因 子 (1 )， 可 得 到 RI-CLPM; 基于 CLPM, 

在 分 解 方 程 中 删除 均值 (px )， 而 保留 截 距 因 子 和 变化 率 因 子 ([1i + (t 一 1)Sx:])， 可 得 到 
LCM-SR; 基于 CLPM， 在 测量 方程 中 保留 测量 误差 (xie )， 在 分 解 方程 中 删除 均值 、 截 距 
因子 和 变化 率 因 子 〈[px + Le + (t 一 1)Sxi])， 在 分 解 方 程 中 保留 恒定 变化 因子 (Gx; )， 可 得 
BJ LCS (HF LCS 是 基于 差 值 的 ， 其 自 回 归 系 数 应 为 Bx = 1 + Bx )。 

此 外 , 有 学 者 证 明了 单 变量 LCM 与 LCS 存在 一 定 的 数理 关联 , 并 提供 了 不 同形 式 下 的 
参数 转换 公式 ， 例 如 ， 对 于 满足 线性 增长 假设 的 单 变 量 LCM, SE (1x;〉 的 均值 和 方差 
等 于 单 变量 恒定 变化 LCS 中 i1 BT IRL TE CC fees ) 的 均值 和 方差 , 其 变化 率 (5%i) 
的 均值 和 方差 则 对 应 LCS 中 的 恒定 变化 因子 (Gxi ) 的 均值 和 方差 , 即 : jx = Ufaan’ Hsx = 


lig, Oh = AME og, = o, (Serang et al., 2019). 


尽管 各 模型 在 结构 上 有 相似 之 处 , 但 不 同 模型 之 间 存 在 明显 差异 , 这 些 差异 有 可 能 导致 
分 析 结 果 及 其 解释 的 角 然 不 同 。 表 1 归纳 了 文中 四 种 交叉 滞后 追踪 模型 的 差异 ， 分 述 如 下 。 
第 一 , 建 模 对 象 不 同 。 CLPM 是 采用 中 心 化 变量 , 即 分 离 出 每 个 时 间 点 的 均值 (在 Mplus 


6 


chinaXiv:202208.00003v1 


中 是 通过 定义 潜 变 量 的 方式 ， 负 荷 都 固定 为 1)，RI-CLPM 是 采用 同时 分 离 出 均值 和 特质 因 


子 后 的 变量 ,LCM-SR 是 采用 去 趋势 的 变量 ， 即 分 离 出 截 距 因 子 和 变化 率 因子 ; LCS 是 基于 
不 同时 间 点 之 间 潜 变量 的 变化 值 来 构建 交叉 滞后 结构 。 


第 二 ， 模 型 预 设 轨迹 不 同 。LCM-SR 主要 用 于 分 析 具 有 线性 或 非 线性 发 展 轨迹 的 变量 


能 提供 关于 变量 发 展 轨迹 特征 的 信息 ; LCS 可 分 析 满 足 指数 型 发 展 轨迹 的 变量 (Mund & 


Nestler, 2019); 其 他 两 个 模型 则 无 法 明确 提供 相关 信息 。 
第 三 , 模型 对 被 试 内 和 被 试 间 效 应 的 分 解 情况 不 同 。 截面 设计 只 能 获得 被 试 间 效 应 的 数 
据 , 而 妃 踪 数据 可 以 同时 反映 被 试 间 差 异 和 被 试 内 过 程 ,区 分 这 两 种 效应 对 于 理解 变量 间 纵 
向 关系 的 本 质 非常 重要 。CLPM 和 LCS 没有 对 被 试 间 差异 成 分 进行 建 模 ， 也 没有 区 分 被 试 
间 差 异 和 被 试 内 过 程 〈Curran et al., 2014; Mund & Nestler, 2019)， 而 RI-CLPM 和 LCM-SR 
则 分 别 通过 分 离 特质 因子 〈 如 je) 和 增长 因子 (如 Jw，Sxi ) 的 影响 ， 区 分 出 了 被 试 间 差 异 。 
如 果 有 足够 理由 认为 在 变量 的 协同 发 展 过 程 中 , 被 试 间 差 异 是 不 可 忽视 的 , 则 研究 者 应 选择 


能 够 区 分 被 试 内 和 被 试 间 效 应 的 模型 ， 从 而 获得 真实 的 变量 间 相互 关系 。 


第 四 ， 模 型 对 于 追踪 数据 的 时 间 点 个 数 要 求 不 同 。Usami A (2019) 指出 ,在 稳定 性 


假设 下 《 即 自 回归 和 交叉 滞后 系数 不 随时 间 变 化 )，CLPM 需要 至 少 两 次 重复 测量 ， 其 他 交 


又 滞后 模型 则 需要 至 少 三 


次 。 


第 五 ， 模 型 适用 的 下 


完 情境 不 同 。 例 如 ， 若 其 中 一 个 变 咎 


量 本 身 并 不 具备 系统 的 变化 趋势 


CUD y 是 阅读 成 绩 , x 是 


可 以 看 出 ,尽管 不 同 模 型 存在 一 定 关 联 , 但 各 个 模型 的 适用 情境 不 同 ， 所 能 提供 的 变量 
发 展 信 息 也 有 所 不 同 (Parker etal., 2015)， 因 此 ， 


个 学 年 中 的 缺勤 天 数 ), 则 不 适合 采用 预 设 发展 轨 迹 的 LCM-SR。 


型 进行 预 设 ， 进 行 模型 比较 和 模型 选择 是 有 必要 的 。 


在 实际 运用 中 ， 如 果 没 有 明确 的 理论 对 模 


表 1 各 交叉 滞后 追踪 模型 的 比较 


建 模 对 象 考虑 测量 ”交叉 滞后 SEH ” 预 设 发 展 时 间 点 
被 试 间 差 试 内 过 适用 | 
(特点 ) 误差 结构 历时 效应 轨迹 oat ANSE 个 数 Wo 
各 时 间 点 潜 变量 自 回归 系 数 、 
CLPM x ~ ~ x =2 \ 存 在 稳定 的 被 试 间 差 异 
(中 心 化 ) 无 en 不 存在 稳定 的 被 试 间 差异 
各 时 间 点 潜 变 量 — 
: 心 化 、 分 离 特质 y " 村 质 dpi > 在 稳定 的 被 试 间 差异 
RI-CLPM (中 心 化 、 分 离 特质 x 无 特质 因子 交叉 滞后 系数 3 存在 稳定 的 被 试 间 差 异 
AT 
各 时 间 点 潜 变量 线性 自 回归 系数 、 
LCM-SR x ~ 以 Pa = 个 变量 均 存 在 系统 的 时 间 趋 
NON gh 长 因子 iesu 3 ”两 人 名 存在 系统 的 时 间 趋 势 
各 时 间 点 之 间 
ALVA AR. 一 个 变量 的 性 水 平 对 于 =A 
LCS irae htt i J 4 y 指数 型 x i 
vies 交叉 滞后 系数 变量 之 后 的 变化 量 的 影响 
(考虑 测量 误差 ) 


注 : 此 处 特 指 双 变量 情况 。 
“VY ”表示 存在 ,“ X ”表示 不 存在 或 无 法 体现 。 
CLPM: 交叉 滞后 面板 模型 ，RILCLPM: 随机 稚 距 交叉 滞后 面板 模型 ，LCM-SR: 结构 化 残 差 潜在 曲线 模型 ，LCS: 潜在 变化 分 数 模型 ， 下 同 。 


4 实证 示例 

下 面 用 一 个 实例 来 说 明 不 同 的 交叉 滞后 追踪 模型 运用 于 同一 批 数据 时 所 得 到 的 分 析 结 
果 的 差异 。 数 据 来 源 于 中 国家 庭 追踪 调查 (China Family Panel Studies, CFPS ) 的 开放 数据 库 。 
拟 探究 主观 幸福 感 和 身体 质量 指数 (Body Mass Index, BMI) 间 的 纵向 关系 ， 主 观 幸 福 感 由 
一 道 问卷 题目 测量 ， 其 要 求 被 试 按照 1 (非常 差 ) BIS GRA) 的 等 级 评价 自己 当前 的 生 
i5: BMI 由 被 试 自我 报告 或 他 人 代 答 的 身高 和 体重 数据 计算 : BMI= 体 重 (千克 ) /身高 〈 厘 
^O 的 平方 。 选 取 三 个 时 间 点 〈T1=2014，T2=2016，T3=2018)， 删 除 异 常 值 并 保留 那些 在 
三 个 时 间 点 都 有 记录 的 被 试 ， 最 终 获 得 来 自 1525 名 被 试 的 有 效 数 据 ， 在 初始 时 间 点 上 被 试 
平均 年 龄 为 20 岁 (范围 为 16-71 岁 )，49% 为 男性 。 

使 用 Mplus 8.0 软件 分 析 《〈 见 附录 )。 此 处 主要 目的 在 于 模型 比较 ， 因 此 均 采 用 默认 变量 
间 存 在 纵向 相互 关系 的 全 模型 ， 且 将 同一 类 型 的 自 回 归 系 数 、 交 叉 滞后 系数 设置 为 不 随时 间 
变化 。 此 外 ， 考 虑 到 LCM-SR 预 设 了 变化 轨迹 ， 先 分 别 对 两 个 变量 进行 单 变量 LCM 分 析 ， 
结果 表明 两 个 变量 都 具备 系统 性 的 时 间 趋 势 : 了 m8=3.80 (p<0.001 ), S 2:0%30e=0.01 (p<0.05 )， 


m 
fiti 


Iz] 


Igwr-20.88 (p<0.001), Sgwur-0.08 (p<0.001)， 适 合用 于 后 续 分 析 。 
不 同 模型 的 拟 合 结果 见 表 2。 当 拟 合 指数 满足 以 下 条 件 : RMSEA 小 于 0.08, SRMR 小 
F 0.08, CFI 大 于 0.90, TLI 大 于 0.90, 可 以 认为 模拟 能 够 较 好 地 拟 合 数据 ( 温 忠 麟 等 , 2004 )。 
可 以 看 出 ，LCM-SR 拟 合 较 差 ， 其 余 模型 的 整体 拟 合 结果 均 较 好 。 
表 2 各 交叉 滞后 追踪 模型 的 拟 合 情 况 


L df RMSEA SRMR CFI TLI 

CLPM 19.68 6 0.04 0.03 0.99 0.99 
LCM-SR 82.79 2 0.16 0.03 0.98 0.83 
RI-CLPM 34.60 6 0.06 0.06 0.99 0.98 
LCS 20.79 7 0.04 0.03 0.99 0.99 


i£: RMSEA-Root-mean-square error of approximation; SRMR=Standard root mean-square residual; CFI= 


Comparative fit index; TLI=Tucker-Lewis index. 


不 同 模型 的 参数 估计 结果 见 表 3。 对 于 同一 批 妃 踪 数 据 ， 各 模型 提供 了 不 同类 型 的 参数 
言 息 。 首 先 ， 对 于 变量 发 展 轨 迹 ，LCM-SR 的 结果 表明 BMI 有 显著 的 线性 增长 趋势 ， 其 截 
距 、 变 化 率 均 显 彰 ， 而 主观 幸福 感 的 变化 率 因子 不 显著 ; 其次， 对 于 自 回归 效应 ，CLPM 和 
RI-CLPM 都 表明 两 个 变量 均 存 在 一 定 的 等 级 稳定 性 , 但 各 模型 所 反映 的 稳定 性 程度 都 不 同 ， 
而 LCM-SR 仪 识别 出 BMI 的 等 级 稳定 性 ，LCS 由 于 模型 较为 复杂 ， 部 分 系数 在 标准 化 后 出 
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现 了 不 恰当 解 (Usami etal., 2019); 最 后 ， 对 于 交叉 滞后 效应 ，CLPM、RI-CLPM fil LCM- 
SR 都 识别 出 了 变量 间 的 历时 性 影响 ， 但 所 得 到 的 作用 方向 不 同 ，CLPM 的 结果 表明 BMI 对 
主观 幸福 感 有 显著 的 历时 性 影响 ， 即 先前 的 BMI 正 向 影响 随后 的 主观 幸福 感 ， 而 RI-CLPM 
的 结果 表明 先前 的 主观 幸福 感 对 随后 的 BMI 有 显著 的 负 向 影响 ，LCM-SR 的 结果 则 表明 先 
前 的 BMI 负 向 影响 随后 的 主观 幸福 感 。 

可 以 看 出 ,在 分 析 同 一 批 数据 时 ， 不 同 模型 所 得 出 的 参数 结果 有 一 定 的 差异 ， 这 是 因为 
不 同 模型 的 构成 不 同 、 前 提 假 设 也 有 所 差异 。 一 方面 ,有 些 模 型 可 能 并 不 适合 分 析 该 批 数据 ， 
例如 LCM-SR 的 模型 拟 合 较 差 ， 说 明 该 模型 可 能 无 法 准确 反映 数据 背后 的 真实 变量 关系 ; 
另 一 方面 ， 由 于 各 模型 的 构成 不 同 ， 参 数 所 反映 的 效应 可 能 是 不 同 的 ， 例 如 CLPM 被 认为 
混淆 了 被 试 间 过 程 和 被 试 内 过 程 ， 因 此 其 参数 结果 可 能 无 法 真实 地 反映 变量 间 的 相互 关系 ， 
而 RI-CLPM 通过 提取 出 稳定 的 被 试 特质 因子 来 分 离 被 试 间 差 异 ， 其 估计 结果 能 够 更 准确 地 
反映 变量 间 的 动态 相互 作用 。 


表 3 各 交叉 滞后 追踪 模型 的 参数 估计 结果 


CLPM RI-CLPM LCM-SR LCS 


I sme = 一 3.85 3.88 一 
NEUES — — —0.02 — 
7BMI a = 20.71" 20.71" 
SBMI 一 — 0.43"* = 
A xen 一 = — 11.55 
ABMI x — = —0.40 
主观 幸福 感 ro >SEMAR om) — 0.65" 0.15" 0.12 x 
主观 幸福 感 rm, 一 主观 幸福 感 r 0.6377 0.18 0.15 x 
BMI cro > BMI» — 0.887 052" 0.68" -0.14 
BMI 2 > BMI m> — 0.9077 0.60 0.61" -0.16 
BMI co 一 主观 幸福 感 m， 0.087 -0.05 -0.37* x 
BMI r 30 5E38/8 e —— 0.08" -0.07 -0.34 x 
主观 幸福 感 o > BMI» -0.01 -0.05* -0.12 0.23 
主观 幸福 感 5 > BMI«a» — -0.01 -0.06™ -0.14 0.30 
了 了 主观 幸福 感 <> [par 一 0.26" = = 


注 : 表 中 路 径 系数 均 为 标准 化 后 的 值 ， 由 于 模型 中 的 等 值 限制 是 对 非 标准 化 路 径 施加 的 ， 因 此 标准 化 路 径 
系数 可 能 会 有 不 相等 的 情况 ; 

LCS 的 路 径 系数 表示 一 个 变量 取 值 对 另 一 个 变量 的 变化 量 的 影响 ; 

LCS 较 复 杂 ， 部 分 路 径 系 数 在 标准 化 后 有 不 恰当 解 〈 用 “x” 表 示 )， 但 整体 模型 拟 合 良好 ， 因 此 仍 报告 
其 他 可 提供 有 用 信息 的 参数 ; 

==, p«0.001; **, p«0.01;*, p<0.05; “<>” 表 示 相 关 . 
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除了 LCM-SR， 其 余 模 型 拟 合 都 较 好 ， 仅 根据 模拟 拟 合 情 况 无 法 确定 最 终 模 型 ， 还 需要 
结合 模型 特点 和 相关 理论 背景 进行 判断。 首先 ， 有 研究 表明 主观 幸福 感 和 BMI 均 存 在 显著 
的 被 试 间 差 异 CBieda etal., 2019; Pereira & Coelho, 2013)， 所 以 排除 CLPM; 此 外 ， LCM- 
SR 拟 合 较 差 ， 说 明 该 模型 不 适合 用 来 反映 所 论 的 变量 关系 ; 最 后 ， 本 研究 的 重点 不 在 于 变 
化 量 ， 且 LCS 过 于 复杂 导致 不 恰当 解 ， 因 此 排除 LCS。 综 上 ， 最终 选择 RI-CLPM 并 得 出 结 
ve: 在 被 试 间 水 平 上 ， 主 观 幸 福 感 与 BMI 存在 显著 正 向 相关 (r=0.26,p<0.05)， 即 主观 幸福 
感 高 的 被 试 群体 倾向 于 有 较 高 的 BMI， 正 所 谓 心 宽 体 胖 ; 在 被 试 内 水 平 上 ， 两 个 变量 均 表 现 
出 一 定 的 等 级 稳定 性 ， 主 观 幸福 感 的 等 级 稳定 性 较 低 (0.15 ~0.18)，BMI 的 等 级 稳定 性 较 证 
(0.52 ~ 0.60); 两 个 变量 间 的 历时 性 影响 : 先前 的 BMI 负 向 影响 随后 的 主观 幸福 感 〈)Bwr 


xs 7-0.05 一 -0.07，P>0.05)， 但 不 显著 ， 先 前 的 主观 幸福 感 负 向 影响 随后 的 BMI Cy ame 
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&-BMI--0.05 ~ -0.06，p<0.05)， 但 影响 微小 。 

5 讨论 

新 兴 追 踪 模 型 多 采用 交叉 滞后 路 径 来 探究 变量 间 的 历时 性 作用 , 因此 各 模型 的 结构 相似 ， 
但 其 中 的 细微 差异 可 能 导致 分 析 结 果 的 不 同 , 尤其 是 在 对 变量 间 影 响 关 系 的 判断 上 可 能 给 出 
在 统计 上 很 不 同 的 结果 ， 甚 至 可 能 得 到 完全 相反 的 方向 ， 因 此 必须 要 选择 适当 的 模型 。 
研究 者 需要 了 解 不 同 追踪 模型 的 适用 情境 ,树立 模型 比较 和 模型 选择 的 意识 ， 并 逐渐 建 


立 追 踪 研 究 的 流程 规范 。 
首先 ,在 研究 设计 阶段 , 要 保证 足够 的 被 试 量 和 时 间 点 , 并 注意 对 测量 
= 一 般 采 用 等 间隔 设计 ， 以 保证 有 意义 的 结果 解释 。 
然后 ， 在 数据 分 析 阶 段 ， 可 根据 数据 特征 (如 时 间 点 个 数 ) 考虑 备 选 模型 ， 多 个 模型 有 
助 于 获得 有 关 变 量 发 展 的 丰富 信息 (Parker et al., 2015). 
最 后 ， 通 过 模型 比较 作出 选择 ,可 从 以 下 几 点 综合 考虑 : (1) 研究 目的 ， 如 果 研 究 者 的 
目的 在 于 描述 变量 的 发 展 轨迹 特征 ， 可 以 选择 LCM 或 LCM-SR， 如 果 目 的 在 于 分 析 变 量 之 
间 的 纵向 相互 作用 机 制 ， 可 以 选择 LCM-SR、CLPM、RI-CLPM 或 LCS; (2) 理论 背景 或 经 
验 , 如 果 有 证 据 表 明 变量 存在 被 试 间 差异 , 需要 从 被 试 间 差异 中 分 离 出 被 试 内 过 程 , 则 LCM- 
SR 或 RI-CLPM 是 较 适 当 的 ; (3) 模型 拟 合 情况 ， 研 究 者 可 借助 各 拟 合 指 数 进行 模型 选择 。 
基于 交叉 滞后 结构 的 追踪 模型 重点 关注 变量 间 的 历时 性 作用 , 可 以 提供 变量 发 展 和 变化 
的 多 种 信息 , 对 于 揭示 变量 间 纵 向 关系 具有 重要 作用 , 能 够 为 变量 间 因 果 关 系 的 论证 提供 证 
据 。 根 据 实际 背景 、 研 究 目 的 和 统计 结果 进行 模型 比较 ， 有 助 于 建立 合适 的 追踪 模型 。 
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本 文 仅 针对 双 变 量 追 踪 模 型 的 全 模型 形式 进行 示例 , 在 实际 分 析 中 , 研究 者 还 可 以 对 全 
模型 进行 删 减 得 到 单 向 模型 或 放宽 稳定 性 假设 的 限制 , 也 可 以 纳入 时 变 或 非 时 变 的 协 变量 做 
进一步 深入 分 析 。 本文 所 述 的 几 种 模型 都 是 在 理论 和 实证 研究 中 得 到 较 多 讨论 和 运用 的 , 此 
外 , 还 存在 其 他 交叉 滞后 模型 ， 如 因子 -交叉 滞后 面板 模型 (Factor-Cross Lagged Panel Model, 


factor-CLPM) (McArdle, 2009) ， 稳 定 特质 自 回 归 特 质 状态 模型 (Stable Trait Autoregressive 


Trait and State Model, STARTS) (Kenny & Zautra, 2001)， 以 及 多 水 平 交 义 灌 后 模型 (Multi- 
level Cross Lagged Model) (Yuetal.,2015) 等 ， 未 来 可 纳入 其 他 交叉 滞后 模型 进行 更 系统 的 


比较 。 
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附录 
Mplus 语句 


Data: FILE = 1.dat; 
Variable: 
names = x1 x2 x3 yl y2 y3; 
Analysis: 
ESTIMATOR=ML; 
ICLPM 
model: 
! define latent variables 
etaxl by x1@1;etax2 by x2@1;etax3 by x3@1; 
etayl by y1(g)l;etay2 by y2(g)l;etay3 by y3@1; 
! autoregressions 
etax2 on etax] (al);etax3 on etax2 (al); 
etay2 on etayl (a2); etay3 on etay2 (a2); 
! crosslagged paths 
etay2 on etax1 (c1);etay3 on etax2 (c1); 
etax2 on etayl (c2);etax3 on etay2 (c2); 
! variance structure 
etaxl;etayl;etaxl with etayl; 
etax2-etax3 (e1);etay2-etay3 (e2); 
etax2 with etay2 (cove);etax3 with etay3 (cove); 
x1-x3;yl-y3; 
! mean structure 
[etax1-etax3 ];[etay1-etay3]; 
[x1-x3(20]; [y1-y3(20]; 


IRI-CLPM 

model: 

ix by x1(21 x2@1 x3@1; 

iy by yI@1 y2@1 y3(21; 

etaxl by x1@1;etax2 by x2@1;etax3 by x3@1; 
etayl by y1(g)l;etay2 by y2@ 1;etay3 by y3@1; 
etax2 on etax] (al);etax3 on etax2 (al); 

etay2 on etayl (a2);etay3 on etay2 (a2); 

etay2 on etax1 (c1);etay3 on etax2 (c1); 

etax2 on etayl (c2);etax3 on etay2 (c2); 
x1-x3@0;y1-y3@0; 

etax 1-etax3;etay 1-etay3; 

Ix; iy; Ix with iy; 

etaxl with ix@0; etayl with ix@0; 

etaxl with iy@0; etayl with iy@0; etax1 with etay1; 


etax2 with etay2 (e1); etax3 with etay3 (e2); 
[x1-x3(20]; [y1-y3(20]; 

[etax1-etax3(20]; [etay1-etay3(2)0]; 

[ix iy]; 


! LCM 

Model: 

il sl | x1@0 x2@1 x3@8 ; 
12 s2 | y1@0 y2@1 y3@2; 
sl with 12; s2 with 11; 


ILCM-SR 

model: 

ix by x1(2)1 x2(2)1 x3@1 ; sx by x1@0 x2(2] x3(22 ; 

iy by y1@1 y2@1 y321 ; sy by y1(20 y2@1 y3(22 ; 

etaxl by x1(2)1; etax2 by x2@1; etax3 by x3@1; 

etayl by y1(g)1; etay2 by y2@1; etay3 by y3@1; 

etax2 on etax] (a1); etax3 on etax2 (al); 

etay2 on etayl (a2); etay3 on etay2 (a2); 

etay2 on etax1 (c1); etay3 on etax2 (c1); etax2 on etayl (c2); etax3 on etay2 (c2); 
x1-x3@0; y1-y3(2)0; etax1-etax3; etayl-etay3; 

etaxl with etayl; 

etax2 with etay2 (e1); etax3 with etay3 (el); 

ix; iy; SX; Sy; 

Ix with sx iy sy; sx with iy sy; iy with sy; 

etaxl with ix@0; etax1 with sx@0; etax1 with iy@0; etax1 with sy@0; 
etayl with ix@0; etayl with sx@0; etayl with ty@0; etayl with sy@0; 
[x1-x3(20]; [y1-y3(20]; 

[etax1-etax3(2)0]; [etay1-etay3(20]; [ix iy]; [sx sy]; 


ILCS 

model: 

etaxl by x1@1; etax2 by x2@1; etax3 by x3(21; 
etayl by yl1@1; etay2 by y2@1; etay3 by y3@!1; 
dx2 by etax2@1; dx3 by etax3@1; dy2 by etay2@1; dy3 by etay3@1; 
gx by dx2-dx3@1; gy by dy2-dy3(2)1; 

etax2 on etax1@1;etax3 on etax2@1; 

etay2 on etayl(g)l;etay3 on etay2@1; 

dx2 on etax1 (p1); dx3 on etax2 (pl); 

dy2 on etayl (p2); dy3 on etay2 (p2); 

dy2 on etax1 (c1); dy3 on etax2 (c1); 

dx2 on etayl (c2); dx3 on etay2 (c2); 

x1-x3 (ux); yl-y3 (uy); 
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etax2-etax3@0;etay2-etay3@0; 
dx2-dx3@0; dy2-dy3@0; 
etax1; etayl; 

g:gy; 

etaxl with etay1; 

gx with etax1; gx with etayl; 
gy with etax1; gy with etayl; 
gx with gy; 

[x1-x3(20]; [y1-y3(20]; 
[etax2-etax3@0]; [etay2-etay3(2)0]; 
[dx2-dx3@0]; [dy2-dy3@0]; 
[etax]l etayl]; 


[gx gy]; 
OUTPUT: STDYX; 


Exploring the longitudinal relations: 


Based on longitudinal models with cross-lagged structure 


Fang Junyan !, Wen Zhonglin ?, Huang Guomin 3 
(! School of Psychology, South China Normal University, Guangzhou, 510631) 
( Center for Studies of Psychological Application / School of Psychology, South China Normal University, Guangzhou, 510631) 


C School of Psychology, South China Normal University, Guangzhou, 510631) 


Abstract A cross-lagged structure usually consists of two kinds of effects, autoregressive effects of the prior level 
of a variable on the current level of itself and cross-lagged effects of the prior level of one variable on the current 
level of another variable. Longitudinal models with the cross-lagged structure are well recognized as powerful 
techniques for revealing longitudinal relations between two variables and laying the foundation of diachronic 
causation. There exist several cross-lagged longitudinal models, while practitioners know little about the association 
and difference among them, which makes it difficult to choose the most proper one. Although these models are 
similar in structure, they may differ in the results of estimation. Thus, it is necessary to get a whole picture of these 
longitudinal models and learn how to compare and choose among them. The present study aims to analyze different 
cross-lagged longitudinal models and compare them, so as to reveal the importance of model comparison and model 
selection and provide strategies to select among models. 

First, we introduce four popular longitudinal models with cross-lagged structure: Cross-Lagged Panel Model 
(CLPM), Random-Intercept Cross-Lagged Panel Model (RI-CLPM), Latent Curve Model with Structured Residuals 
(LCM-SR), and Latent Change Score Model (LCS). Then, we clarify the similarities and associations among them. 
Next, we discuss their differences in various aspects. Finally, we conduct an empirical study to illustrate the 
procedure of model selection. 

Results show that: (1) these models are very similar in the model configuration because they all analyze 
diachronic relations by the cross-lagged structure; (2) CLPM can transform into RI-CLPM, LCM-SR and LCS under 
certain conditions; (3) different models focus on different developmental characteristics and each of them can 
provide valuable information on the change process; (4) these models could give different estimation results when 
applied to the same data set, which may induce different conclusions. 

We summarize several reference points for selecting a proper longitudinal model in practice: (1) research 


purpose. If researchers are interested in characterizing the development trajectories, then LCM or LCM-SR is 


preferred; (2) theoretical knowledge and empirical experience. If there is sufficient evidence showing that the within- 
person process should be separated from between-person difference, then LCM-SR and RI-CLPM could be 
considered; (3) the model fitting. Several model fit indices can be used. 

In summary, longitudinal models with cross-lagged structure play an important role in revealing longitudinal 
relations between psychological constructs. These models are similar in configuration but vary in modeling basis, 
premises and data requirements, which may give rise to distinct estimation results and conclusions. Researchers 
should understand the association and differences among them with considerable insight into model comparison and 
model selection. It is advisable to try different reasonable models and choose the most proper one for the exploration 
of longitudinal relations. 


Keywords longitudinal relation, cross-lagged, longitudinal model, model selection 


